arXiv:1504.03835vl [cond-mat.quant-gas] 15 Apr 2015 


Proposed method to realize the p-wave superfluid state using a 
s-wave superfluid Fermi gas with a synthetic spin-orbit interaction 

T. Yamaguchi and Y. Ohashi 
Department of Physics, Keio University, 

3-14-1 Hiyoshi, Yokohama 223-8522, Japan 
(Dated: April 16, 2015) 

Abstract 

We theoretically propose an idea to reach the p -wave superfhiid phase in an ultracold Fermi 
gas. The key of our idea is that the pairing symmetry of a Fermi superfluid is fully dominated 
by the symmetry of the superfluid order parameter, which is essentially given by the product of a 
pair amplitude and a pairing interaction. Noting this, in our proposal, we first prepare a p -wave 
pair amplitude by, not using a p -wave interaction, but using the phenomenon that a p-wave pair 
amplitude is induced in an s-wave superfluid Fermi gas with antisymmetric spin-orbit interaction. 
In this case, although the system is still in the s-wave superfluid state with the s-wave superfluid 
order parameter, when one suddenly replaces the s-wave interaction by an appropriate p -wave one 
(which is possible in cold Fermi gases by using a Feshbach resonance technique), the product of 
the p -wave interaction and the p -wave pair amplitude that has already been prepared in the spin- 
orbit coupled s-wave superfluid state immediately gives a finite p -wave superfluid order parameter. 
Thus, at least just after this manipulation, the system is in the p -wave superfluid state, being 
characterized by the artificially produced p -wave superfluid parameter. In this paper, to assess our 
idea, we evaluate the p -wave pair amplitude in a spin-orbit coupled s-wave superfluid Fermi gas 
at T = 0. We determine the region where a large p -wave pair amplitude is obtained in the phase 
diagram with respect to the strengths of the s-wave pairing interaction and the spin-orbit coupling. 
We also discuss the accessibility of this optimal region on the viewpoint of the superfluid phase 
transition temperature. Since the achievement of a p -wave superfluid Fermi gas is one of the most 
crucial issues in cold atom physics, our proposal would be useful for this exciting challenge. 

PACS numbers: 03.75.Ss, 03.75.-b, 67.85.Lm 
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I. INTRODUCTION 


:Q and 6 Li {2 4 


Since the achievement of the s-wave superfluid phase transition in 4 K 
Fermi gases, the possibility of the p-wave superfluid state has extensively been discussed 


as 


the next challenge in cold Fermi gas physics fs-IlQl]. Although no one has succeedec 


in this attempt, p- wave Feshbach resonances have already been discovered in 10 K 


20 


23 


and 6 Li 


case 


29 


24 


26J. Thus, once a p-wave superfluid Fermi gas is obtained, as in the s-wave 


, one can systematically study this unconventional Fermi superfluid from the 


weak-coupling regime to the strong-coupling limit, by adjusting the threshold energy of 


a p- wave Feshbach resonance 27^ 


28|. Since this pairing state has also been discussed in 


various fields, such as liquid 3 He[40H42]. metallic superconductivity |43h45] . as well as neutron 
star 46| , the realization of a tunable p-wave superfluid Fermi gas would also give a big impact 


in these research fields. 

In the current stage of research toward the realization of a p-wave superfluid Fermi gas, 
one serious difficulty is that, although a p-wave attractive interaction is necessary for this 


pairing state, it also causes the so-called three-body loss 


relaxation 


Bid 


23], leading to 


very s. 


rort lifetime of p-wave pairs 


48], as well as the dipolar 


49]. Thus, although one can 


prepare p-wave molecules(23. I50 h 52] . they are soon destroyed, before the p-wave condensate 


grows. Thus, at this stage, it is crucial to overcome this difficulty. 

The purpose of this paper is to propose a possible route to reach the p-wave superfluid 
phase in an ultracold Fermi gas. To explain our idea in a simple manner, we first note 
that any Fermi superfluid is characterized by a superfluid order parameter A aa /(p), being 
consisting of the product of the pairing interaction U(p, p') and the pair amplitude ( p ) = 

(c Pi(T C-p,o-') as 


^a,a'{jP) ^ ^ U{p, p ) (Cp/ !(r C_p/ j( j/) 

P' 

= Ul P ^ ^ Tp'^V.ct' (p ) • (1) 

p' 

Here, c Pi(T is the annihilation operator of a fermion with momentum p and (pseudo) spin 
a = j\ i- In the second line in Eq. ([T]) , we have assumed the separable interaction U(p, p') = 
U 7 p 7 p/ (where U is a coupling constant and 7 P is a basis function). When the pair amplitude 
< h( 7 ,( 7 , (p) has the same symmetry as that of the basis function y p , we obtain a finite value of 
the superfluid order parameter A cr 0 -/(p) oc 7 p . 
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In a Fermi superfluid, the pair amplitude < f > 0 -, 0 -'(p) is usually produced by the pairing 
interaction U(p,p') of the system. Thus, current experiments usually deal with a Fermi gas 
with a p- wave pairing interaction from the beginning. Although this approach seems the 
shortest way to reach the p- wave superfluid state, at present, no one has succeeded in this 
attempt, because of the above mentioned serious problem. 

Instead of this ordinary approach, we consider an alternative route. That is, we first 
prepare only a p-wave pair amplitude in a system with no p- wave interaction. Of course, even 
when we prepare this quantity, the system is still not a p- wave Fermi superfluid, because 
the p-wave superfluid order parameter is absent due to the vanishing p-wave interaction. 
However, during this preparation, various difficulties originating from the p-wave interaction 
can be avoided. After this preparation, when the s-wave pairing interaction is suddenly 
replaced by a p-wave one U p (p,p') = Up'jp'yp' (where the p-wave basis function y p is chosen 
so that the momentum summation in Eq. (jT]) can be finite), the product of this interaction 
and the p-wave pair amplitude that has been prepared in advance immediately gives a finite 
value of the p-wave superfluid order parameter A a ,o'(p) oc 7 P . Thus, at least just after this 
manipulation, the system is regarded as a p-wave Fermi superfluid, by definition. 

In a sense, the superfluid order parameter is artificially produced in our approach, so 
that the resulting p-wave superfluid Fermi gas would initially be in the non-equilibrium 
state. In addition, when the p-wave interaction is turned on, the above mentioned problem 
occurs, so that the p-wave superfluid phase might eventually be destroyed. However, the 
p-wave superfluid order parameter would continue to exist for a while, during the decay 
and relaxation of the system. This is quite different from the ordinary approach, where the 
system is destroyed before the reach of the p-wave superfluid phase. 

To prepare a p-wave pair amplitude in our approach, the recent artificial gauge field 


would be useful 


53 


61]. This sophisticated technique enables us to introduce an antisym¬ 


metric spin-orbit interaction to an ultracold atom gas, leading to the broken spatial inversion 


symmetry. In this case, the parity is no longer a gooc 


symmetry of the Fermi superfluid order parameter 62 


quantity to classify the spatial pairing 


63]. Since the pair wavefunction must 


be antisymmetric with respect to the exchange of two fermions, the parity mixing naturally 


leads to the admixture of the spin singlet and spin triplet states|62j, |63|. As a result, even in 


a s-wave Fermi superfluid, the pair amplitude may have the spin-triplet component. Thus, 
using this, one may prepare a p- wave pair amplitude without relying on a p- wave interaction. 
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Indeed, Refs. [64], [65] explicitly show that the p -wave pair amplitude is induced in a s-wave 


superfluid Fermi gas in the presence of the Rashba spin-orbit interaction 
the Hamiltonian, 


661 ]. described by 


H, 


Rashba 


= A 


Rashba 




v a a,a ']c^ c / 
y i J y u x J i 


( 2 ) 


where <7j(i = x,y,z) are Pauli matrices acting on spin space and AR as hba is a spin-orbit 
coupling constant. At present, a simple antisymmetric spin-orbit interaction has been re¬ 


alized in 6 Li 


and 40 K 


59] Fermi gases anc 


spin-orbit couplings have also been proposed 


various ideas to synthesize more complicated 


60 


61]. 


Once a p-wave pair amplitude is induced in a s-wave superfluid Fermi gas, the p-wave 
superfluid order parameter can be produced by the sudden replacement of the s-wave in¬ 
teraction with a p-wave one, by tuning an external magnetic field from a s-wave Feshbach 
resonance field to a p-wave Feshbach resonance field. After this manipulation, the s-wave 
superfluid order parameter immediately vanishes, because of the vanishing s-wave interac¬ 
tion. Although the s-wave pair amplitude still remains, since the pairing symmetry of a 
Fermi superfluid is determined by the symmetry of the order parameter, this system is in 
the p-wave superfluid state. 

We note that the admixture of the spin-singlet pairing and spin-triplet pairing by 
an antisymmetric spin-orbit interaction has also been discussed in non-centrosymmetric 


superconductivity 


63]. For example, Ref.[67| points out the importance of this admix¬ 


ture to understand the anomalous behavior of the penetration depth observed in the non- 
centrosymmetric superconductor L^PtE^. 

Our idea is somehow related to the proximity effect in a superconductor-normal metal 
(S-N) junction 68]. In this case, the pair amplitude in the S side penetrates into the N side, 
so that, when it couples with an interaction existing in the N side, a finite superconducting 
order parameter appears in the N side. Since the pair amplitude in the N side is fully supplied 
from the S side, even when the interaction in the N side is repulsive , this proximity-induced 
superconducting state is obtained 69]. 

In this paper, to assess our idea, we investigate a s-wave superfluid Fermi gas with an anti 


symmetric spin-orbit interaction. Using the strong-coupling theory developed by Eagles 
and Leggett 


30 


29] 


31] (which is sometimes referred to the BCS-Leggett theory in the litera¬ 


ture), we examine how large a p -wave pair amplitude is induced by a spin-orbit interaction 
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at T — 0. We determine the region where a large p -wave amplitude is obtained, in the phase 
diagram with respect to the strengths of the s-wave interaction and the spin-orbit coupling. 
We also examine the accessibility of this “optimal region” within the current experimental 
technique. For this purpose, employing the strong-coupling theory developed by Nozieres 
and Schmitt-Rink (NSR) 32], we calculate the superfluid phase transition temperature T c 
around this region. We briefly note that T c in a spin-orbit coupled ultracold Fermi gas has 
recently been discussed in the BCS (Bardeen-Cooper-Schrieffer)-BEC_(Bose-Einstein con¬ 


densation) region within the framework of a T-matrix approximation 


fluctuations are treated within the static approximation 


30 , 


where pairing 


3d ]. In this paper, we do not em¬ 


ploy the static approximation, but fully take into account dynamical properties of pairing 
fluctuations within the NSR scheme. 

For the time evolution of the system after the introduction of a p-wave interaction, we 
need to deal with a non-equilibrium Fermi gas, which we will separately discuss in our future 
paper. 

This paper is organized as follows. In Sec. 11, we explain the outline of our formulations 
at T — 0, as well as at T c , in the BCS-BEC crossover region. Here, we also introduce the 
condensate fraction, as a useful quantity to evaluate the magnitude of the pair amplitude. 
In Sec. Ill, we discuss how large a p -wave pair amplitude is induced at T = 0, when the 
inversion symmetry of the system is broken by a spin-orbit interaction. We clarify the 
condition to obtain a large p -wave pair amplitude. We also calculate the superfluid phase 
transition temperature T c , to examine the accessibility of this condition within the current 
experimental technique. Throughout this paper, we set h = k-& = 1, and the system volume 
is taken to be unity, for simplicity. 

II. FORMULATION 


We consider a spin-orbit coupled uniform two-component Fermi gas with an s-wave pair¬ 


ing interaction. 
Leggett theory 


29 


'o explain our formulation at T = 0 within the framework of the BCS- 


31], as well as the formulation at T c within the framework of the NSR 


theory[321]. in a unified manner, the functional integral formalism is convenient 33j, [72]. The 
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partition function in this formalism is given by 


Z = 




-5 


(3) 


where the action S has the form, 


S= dx\j2^ 


(x 


dp 2 

dr +2 m ^ . 


a, a' 

(4) 


Here, x = (r, r) and f dx — dr f dr , where fd = 1/T. 'h (T (r,r) and d r 0 -(r,r) are a 
Grassmann variable and its conjugate, describing Fermi atoms with the atomic mass m and 
pseudospin a =t, I) specifying two atomic hyperfine states, p = — i V is the momentum 
operator in real space and /i is the Fermi chemical potential. The s-wave pairing interaction 
— U s (< 0) is assumed to be tunable by an s-wave Feshbach resonance. As usual, we measure 
the interaction strength in terms of the s-wave scattering length a s , given by 

47ra, U q 


m 


1 - u s EJ 


Pc _J_ ’ 
2e v 


(5) 


where £ p = p 2 /(2m) and p c is a cutoff momentum. 


The antisymmetric spin-orbit interaction h so = {h^ } with linear-momentum depen¬ 


dence in Eq. (j3J) generally has the form 73476] 


hso ^ ^ Pi A,;jCr ? , 
hi 


( 6 ) 


where the 3 x 3-matrix A = {Ajj} ( i,j = x,y,z ) describes spin-orbit coupling strengths. 
However, this paper deals with the simpler version, A = diag[Aj_, A_l, A-], that is, 


hso A _ \j) x d x T Pydy\ T X z Pzh z 


(7) 


Here, we take Aj_,A 2 > 0 without loss of generality. Although Eq. (J7|) cannot describe 
all possible spin-orbit interactions, it still covers some typical cases that have recently been 
discussed in cold atom physics. The Rashba-type spin-orbit interaction in Eq. (J2]) is obtained 
by setting A z = 0 and rotating the momentum space by 7 t/2 around the p z axis. The single¬ 
component spin-orbit interaction h so ~ p x a yi which has recently been synthesized in 6 Li[58| 
and 40 K 59] Fermi gases, is also obtained by setting A x = 0 and rotating the momentum 
space by 7 t/ 2 around the p y axis, which is followed by the rotation in the (<t x , a y , a z ) space 
by 7 t/ 2 around the d x axis. 
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As usual, we introduce the Cooper pair field A(x) and its conjugate A*(a;), using the 


Hubbard-Stratonovich transformation 
to v h (T (x) and fy a (x), one has 


33]. Carrying out the functional integrals with respect 


Z = I VA*VAe~ SesAA *\ 


( 8 ) 


where the action S' e g(A, A*) is given by 


S eS (A,A*) = f dx 1 -^^- - ^Tr ln[—G _1 ] 


(9) 


Here, 


G 1 (x,x') = 


( _l_ 

dr 


Pi 

2 771 


—ia y A*(x) 


ha 


_d_ 


ia y A(x) 
2 


P- 


-p 


- h* 

1 t'o, 
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5(x — x') (10) 


V ~~ y v ' dr ' L2 m 

is the inverse of the 4 x 4-matrix single-particle thermal Green’s function G(x,x') = 
— (T r {'h(x) v hi(x , )}) in the operator formalism, where 

^ v h t (x) ^ 


4/(x) = 


is the four component Narnbu field 


(77, 78]. 


'hj(x) 

V *!(*)7 


( 11 ) 


A. p-wave pair amplitude in a spin-orbit coupled s-wave superfluid state at T = 0 


We eva 
Leggett 29 


uate the p -wave pair amplitude at T = 0 within the framework of the BCS- 


31] theory. In the functional integral formalism, this strong-coupling theory 


is simply reproduced by approximately evaluating the functional integral in the partition 
function Z in Eq. (]8]) by the value at the saddle point solution (A s ), which is determined 
from the equation, 


0 = 


<5A*(x) 7 A(x)=A*(x)=A s U s 4/5 


+ Jr Tr 1A/ + *pJ'jjG MF (p, iw„) 


( 12 ) 


P,UJ n 


where c o n is the fermion Matsubara frequency and pj (j = x, y, z) are Pauli matrices acting on 
particle-hole space. (Note that &j ( j = x,y,z ) act on spin space.) In Eq. (fl2]h G Mt (p, iix n ) 
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is the Fourier transformed Green’s function in the mean-field BCS approximation, given by 


G MF (p,iu n ) = - 


iu n ~ [£ P + P\ ■ r}p z - p y cr y A s 

1 ^ iu n + [f p + p\ ■ r]p z + A s p y dy 


A£ 


a==b 


ui + (j E: 


1 To; 


Pa ■ r 

IpaI 


(13) 


where £, p = e p — p = p 2 /(2m) — p is the kinetic energy of a Fermi atom, measured from the 
Fermi chemical potential p. E™ = yj (ip) 2 + A^ describes the Bogoliubov single particle exci¬ 
tations, where = i P +a\p x \ with p x = (A ± p x , A ± p y , A z p z ). In Eq. (H3 If = (p z d x , & yi p z a z ) 
is the spin operator under the four component Nambu representation 77,(78]. Summing up 
the Matsubara frequencies in Eq. (fT2lh we obtain the BCS gap equation in the presence of 
spin-orbit interaction. At T = 0, it has the form, 

U„ x— 1 


1 = 


£ 


2 ^ 2FA 

p,a=± v 


(14) 


Eliminating the ultraviolet divergence from Eq. (TT4|) . we obtain 30. ]33] 


1 = 


47ra« 


m 


Z^L2^2E“ 2e„. 

a=± V 


(15) 


We solve the renormalized gap equation (fT5]h together with the equation for the number 
N of Fermi atoms, which is obtained from the mean-field thermodynamic potential = 
TS eS ( A s , A s ) as 


1 y 

dp 2/3 ^ 


p z G Mb (p, iuj n )e lpz 


£ 


i--£ — 

Q! = ± ** . 


(16) 


P,OJn P 

(where 8 is an infinitesimally small positive number), to self-consistently determine A s and 

p. 

The pair amplitude 4>(p, S, S z ) with total spin S and its ^-component S z is obtained from 
the off-diagonal components of the Green’s function G MF (p, ioj n ) in Eq. (fT3|) . Noting that 
G MF (p,r) in the operator formalism has the form 


G m *( P ,t) = -( T r { 


( c p,tO) \ 
c P,t( r ) 
C -P,t( r ) 

V c -p,4.( r ) / 


( c p,f(0), c p,4_(0), c -p,t(0), c -p,t(0))}), 


(17) 
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one finds, 


$0,0,0) = - 


(c P)t C- p4 ) - (c p4 c_ Pit ) =~ [Gf/(p, r = 0) - G^ f O, t = 0) 


1 

‘2^2Eg’ 

a=± * 


$0, !, !) = (cp,t c - P ,t) = ~ G iz O, t = 0) = 

$0, i,o) 


A_l Px - ip y aA 

2 Ipa| 


E 


2 A"’ 


a=± P 


2 L 


( c p ,t c - p4 ) T ( c p4 c - P ,t) 


2 L 


G*14 O, T = 0) + G 23 ( P , r = 0) 

A z p z «A S 


$0, 1 , - 1 ) = (c p4 c_ p4 ) = -G 24 (p, t = 0) = 


2 OaI^ 2 ^’ 

A_L Px + ip y «' 


Pa 


E 


=± 2E p 


(18) 


Equation ([TJ) clearly shows that the antisymmetric spin-orbit interaction induces the p-wave 
pair amplitudes <E>(p, S = 1, S z = ±1, 0). However, we emphasize that the p-wave superfluid 
order parameter is still absent due to the vanishing p -wave pairing interaction. The system 
is thus in the s-wave superfluid state which is characterized by the s-wave order parameter, 


At t 


A s = ~U S ^ 0.0). 


(19) 


P 


0 , we suddenly replace the s-wave pairing interaction by the p -wave one 



H p -wave Up P ' P C v+a/2.a C ^~T>+a/2.a' C -p'+<l/‘ 2 ,<T lC P'+<l/‘ 2 ,<r- ( 2 0) 

P,P',Q 

Then, while the s-wave order parameter immediately vanishes due to the vanishing s-wave 
interaction, a p-wave order parameter A p (p, S z ,t > 0 ) become finite, so that, by definition, 
the system is in the p-wave superfluid state. In particular, the p-wave order parameter just 
after this manipulation (t = + 0 ) is simply given by the momentum summation of the product 
of the p-wave pairing interaction U p and a p-wave pair amplitude $(p, S = 1, S z )\ t =s that 
has already been prepared in the s-wave superfluid state. For example, when = 0 and 
A 2 7 ^ 0, one has 


A p (p,S z = 0,t = + 0 ) = -U p p z ^2p' z ^(p’, l,0)|i =_ 5 oc p z . ( 21 ) 

p 1 


We note that, although all the pair amplitudes &(p,S,S z ) in Eq. (jT 8 l) are proportional 
to the s-wave order parameter A s , it does not mean that they immediately disappear, when 
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A s vanishes. As well known in the proximity effect 68, [69], even when the pair amplitude 
penetrates into the normal metal with no interaction, it continues to exist until they are 
destroyed by external perturbations, such as thermal fluctuations and impurity scatterings. 
Thus, in the present case, the pair amplitudes in Eq. (ITS]) should be regarded as the “initial 
values”, in considering the time evolutions of <h(;p, S, S z ,t > 0) after the sudden change of 
the pairing interaction from the s-wave one to the p-wave one. 


To evaluate the strength of the p-wave pair amplitude, the condensate fraction 


30 


791-482] 


is a useful quantity. It is deeply related to the pair amplitude and physically describes the 
number of Bose-condensed Cooper pairs. In addition it has widely been used in detecting the 


superfluid phase transition in ultracold Fermi ga ses [1 
case, the total condensate fraction Ah is given by { 31 ] 


-4], In the present spin-orbit coupled 


ivj = - Y, (Cj.,I 2 = N C (S = 0 ,S Z = 0) + Y, V(S = 1 ,SJ. 


( 22 ) 


p,(T,cr' 


s z =-1 


Here, 


K(S = O.S. = 0) = £ |f(p,0,0)| 2 = ^E(E A"' 2 


16 ^ E‘ 

P OL=± V 


(23) 


is the s-wave condensate fraction. Equation 
N C (S = 1, S z = ±1, 0), given by 


also involves the p-wave components 


iv„(M) = l5H«( P .i,i)i 2 = §^ Ai(p ' + h ) 


JVc(i.o) = J>(p,i,o)| 2 = 4f£ A ' ri 


16 ^ |pa | 2 


IpaI 2 

E 


S' ’ 


a 


E° 

a=± P . 


%(1,-D = 


Substituting Eqs. 


and 


into Eq. 


32 

p 

, one has 


Pa 


v Of 

E° 

=± p 


(24) 


A 2 


8 


p,a=z t 


1 

E° 


(25) 


In Sec. Ill, we will numerically calculate N C (S — 1 ,S Z — ±1,0) to examine how large the 
p-wave components are induced by a spin-orbit interaction. 
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B. Superfluid phase transition temperature and effects of spin-orbit interaction 


To evaluate T c and effects of a spin-orbit interaction in the BCS-BEC crossover re¬ 
gion, we include pairing fluctuations within the framework of the strong-coupling theory 
developed by Nozieres and Schmitt-Rink 32], In the functional integral formalism, the 
T c equation in the NSR scheme is immediately obtained from the saddle point condition, 
5S e ff/5A*(x) |a(x)=a*(z)=o = 0 33]. After the renormalization, it is given by, 


1 = 


47ra„ 


m 


VU V-i-tanh^-- — 
| 2 2T“ 2T 9^ 

Q = ± ’P 


2e p i 


(26) 


As usual, we solve the T c equation f[26|) . together with the equation for the number N of 
Fermi atoms, to self-consistently determine T c and the Fermi chemical potential /i. The NSR 
number equation includes pairing fluctuations within the Gaussian fluctuation level, which 
is derived from the identity N = — <90nsr/< 9/^ The NSR thermodynamic potential f^sR is 
obtained by expanding the action S e g- in Eq. (|9]) around A(x) — 0 to the quadratic order, 
which is followed by carrying out functional integrals with respect to A(x) and A*(a;). The 
result is 


d 


N = N free ~T—J2 In 


1 + 


47ra, r 


m 


If {q,iv n 


E — 

-^2e„ 




(27) 


q,Vn p 

where we have eliminated the ultraviolet divergence by employing the renormalization 
prescription 33] and 6 is infinitely small positive real number. In Eq. (127)1 . Nf ree = 
Spa=±/(^p) is ^ ie number of Fermi atoms in the absence of the pairing interaction U s , 
where f(x) is the Fermi distribution function. The second term in Eq. (1271) describes effects 
of pairing fluctuations, where 


n (q,iu n ) 


1 1 f(£p+q/ 2) f(£p-q/ 2) r ,(p\ + q\/ 2 ) ■ (px — Qa/ 2 )" 

%*!* +&*/*-*» [ a |PA + W 2 ||PA-W 2 | j 


(28) 


is the lowest-order pair-correlation function in terms of the pairing interaction (where v n is 
the boson Matsubara frequency). 
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FIG. 1: (Color online) Calculated intensity of pair amplitude §(p, S, S z ) at T = 0. We set Ax = 0 
and \ z /vf = 1 (where vp = kp/m is the Fermi velocity with kp being the Fermi momentum of a 
free Fermi gas in the absence of spin-orbit interaction). The upper and lower panels show the p- 
wave component <b(p, S = 1 ,S Z = 0) and the s-wave component <b(p, S = 0 ,S Z = 0), respectively, 
(al) and (bl) ( kpa s ) _1 = —1. (a2) and (b2) ( kpa s ) _1 = 0. (a3) and (b3) ( kpa s ) _1 = 1. In these 
panels, we take p y = 0. 


III. p-WAVE PAIR AMPLITUDE INDUCED BY ANTISYMMETRIC SPIN- 
ORBIT INTERACTION 


A. single-component spin-orbit interaction (Ax = 0, A 2 / 0) 


Figure ED shows the pair amplitude $(p,S,S z ) in a spin-orbit coupled s-wave superfluid 
Fermi gas at T = 0 in the case of Ax = 0 and A 2 ^ 0. Besides the s-wave pair amplitude 
shown in Figs. [U(bl)-(b3), Figs. [Q(al)-(a3) show that this single-component spin-orbit 
interaction induces the p -wave component with (S', S z ) = (1, 0). While the s-wave component 
does not change its sign (See panels (bl)-(b3).), panels (al)-(a3) clearly show the p z -wave 
symmetry, as expected from the expression for <h(p, S — 1, S z — 0) in Eq. (1T81) . 

In the weak-coupling BCS regime, one sees in Figs. [D(al) and (bl) that both the p -wave 
and s-wave pair amplitudes are large around two circles. This is simply because, in the 
absence of the s-wave pairing interaction U s , the present spin-orbit interaction gives two 
single-particle dispersions 


_ Pi , (Pz±m \ z ) 2 _ _ 

p 2m 2m ^ Z1 


(29) 
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where p\ = p 2 x + p 2 and p z — ji + m\ 2 /2 is an effective Fermi chemical potential. These 
bands give two Fermi surfaces that are centered at p — (0, 0, ±m\ z ) with the radius being 
equal to the Fermi momentum kp in the absence of the spin-orbit interaction. In the weak- 
coupling BCS regime, since atoms near these two Fermi surfaces dominantly contribute 
to the Cooper-pair formation, the pair amplitude becomes large around them, as seen in 
Figs. OJal) and (bl). Since atoms away from the Fermi surfaces also contribute to the pair 
formation when the pairing interaction is strong, the circular structure is obscure in Figs. 
QJa2) and (b2). In the strong-coupling BEC regime, the effective Fermi chemical potential 
p z is negative (See Fig. [21(c).), so that the Fermi surface no longer exists. As a result, the 
circular structure disappears in Figs. [U(a3) and (b3). 

To evaluate the p -wave pair amplitude in a quantitative manner, we conveniently consider 
the condensate fraction N C (S, S z ) at T = 0. In the BCS side ((/q?as ) -1 ^ 0), Fig. [21(a) shows 
that the p -wave component N C (S = 1, S z = 0) increases with increasing the interaction 
strength, reflecting the increase of the magnitude of the s-wave superfluid order parameter A s 
shown in Fig. [21(b). However, N C (S — 1, S z — 0) decreases in the BEC side ((kpds)- 1 > 0), 
although the s-wave condensate fraction N C (S — 0, S z — 0) continues to increase. In this 
regime, the coupled equations (1T51) and (1X6]) give, 


so that one obtains 


A, 

p 
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3Tr(kpa s ) F ’ 

- 7 - mX 2 , 

2 ma 2 2 25 


iYc(1 ’ 0) = (M - )2 ’ 
N 

Ac(0,0) = —— AT C (1, 0). 


(30) 


(31) 


This means that the strong-coupling BEC limit ((/cpa s ) _1 —> 00 ) may be simply viewed as 
a gas of N/2 s-wave molecules, even in the presence of spin-orbit interaction. 

When the spin-orbit interaction becomes strong, the magnitude of p-wave component 
N C (S — 1, S z — 0) increases, as shown in Fig. [3](a). However, in the present case of single¬ 
component spin-orbit interaction, we should note that the total condensate fraction IV*, as 
well as the s-wave superfluid order parameter A s , do not depend on the spin-orbit coupling 
strength A,. Indeed, when A^ = 0 and A, 7 ^ 0 , the coupled equations (fl5l) and (TT 6 |) are 
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FIG. 2: (Color online) (a) Calculated condensate fraction N C (S,S Z ) in a spin-orbit coupled s-wave 
superfluid Fermi gas at T = 0. We take \ z /vy = 1 and Aj_ = 0. In this single component case, one 
finds that iV c (l,±l) = 0. N* is the total condensate fraction in Eq. ([22]) . (b) s-wave superfluid 
order parameter A s , normalized by the Fermi energy ep = fcp/(2m). (c) Effective Fermi chemical 
potential fi z = n + m\ 2 z /2. 


reduced to the ordinary BCS-Leggett coupled equations for a s-wave Fermi superfluid with 
no spin-orbit interaction , as 


1 


47ra s 

m 


£ 




N 



(32) 


where £ p = £ p — fi z . Thus, the self-consistent solutions (A s , /x~) shown in Figs. [2(b) and 
(c) are independent of X z . As a result, the total condensate fraction iV* in Eq. (|25|) is also 
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FIG. 3: (Color online) (a) Calculated p -wave condensate fraction N C (S = 1, S z = 0) at T = 0 in the 
case of single-component spin-orbit interaction (A_l = 0). (b) Peak position of N C (S = 1 ,S Z = 0) 
evaluated from the result shown in panel (a). The background intensity shows T c . (c) Calculated 
T c . The inset shows the effective Fermi chemical potential p z = p + m\ 2 /2. 


A,-independent as 


± ' C 

N 


E 


37tA, 


A? 


W V % + idV2ef- 


\j Hz + 


(33) 


where £ F is the Fermi energy for a free Fermi gas without spin-orbit interaction. As a result, 
for a given interaction strength (fc F a s ) _1 , the s-wave condensate fraction N C (S — 0, S z — 
0) = iV* — N C (S = 1, S z = 0) becomes small when the spin-orbit interaction becomes strong, 
in spite of the fact that the s-wave superfluid order parameter A s remains unchanged. 
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Evaluating the peak position of the p-wave condensate fraction N C (S — 1, S z — 0) from 
Fig. 0(a), we obtain Fig. 0(b). Recently, the single-component spin-orbit interaction with 
0.5 < \ z /v F < 1 has been realized in a 40 K Fermi gas[59]. Keeping this in mind, we find 
from Fig. 0(b) that the region 


( k F a s ) ~ 0, \ z /v F ~ 1, 


(34) 


is suitable for the preparation of large p- wave pair amplitude. As an example, at ( k F a s ) 1 = 0 
and X/v F = 1, one obtains 


N C (S = 1,S Z = 0) 
N 


~ 0.07. 


(35) 


In this case, just after the sudden change of the pairing interaction from the s-wave one to 
the p- wave one in Eq. (12UD . we expect the p- wave superfluid state with the p- wave condensate 

;e with A p (p,S i = 0) oc p z is 


fraction being equal to Eq. (I35]h This p- wave superfluid sta 
just the so-called polar phase discussed in superfluid 3 He 
summarize this p- wave superfluid order parameter as 


41 


42], As usual 41443]. one can 


A(p) = 


(36) 


A p (p,S z = l) A p (p,S z = 0) \ I 0 Pz 

A P (p, S z = 0) Ap(p, S z = -1) J y p z 0 

where we have ignored the unimportant factor that is nothing to do with the pairing sym¬ 
metry in the last expression. Since the polar state has not been realized in liquid 3 He, a 
single-component spin-orbit coupled s-wave superfluid Fermi gas would be useful for the 
realization of this unconventional pairing state. 

To actually prepare the p- wave pair amplitude in the parameter region in Eq. (I34|h the 
s-wave superfluid phase in this regime must be experimentally accessible. In this regard, we 
note that T c is also A^-independent in the present case. Indeed, when A_l = 0 and A 2 ^ 0, 
we can completely eliminate the A 2 -dependence from the T c -equation (1261) as 


1 = 


4,r “‘WpLtanhi^- — 

| oe 9 T 2e r 


m 


2£ P 


(37) 


where is given below Eq. (J32J). In the same manner, the non-interacting part Af ree , as well 
as the pair correlation function Il(q, fz/ n ), in the number equation (1271) can be also written 
in the A 2 -independent forms as, respectively, 


A^free 


AZ). 


(38) 
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FIG. 4: (Color online) Intensity of the pair amplitude <h(p, S, S z ) at T = 0 in the case of Ax/uf = 1 
and A- = 0. The upper and middle panels show Re[<h(p, S = 1, S z = 1)] and Im[4>(p, S = 1, S z = 
1)], respectively. The lower panels show the s-wave pair amplitude 4>(p, S = 0, S z = 0). (al)-(cl) 
(/cfOs) - 1 = — 1. (a2)-(c2) = 0. (a3)-(c3) (&Fa s ) _1 = 1. In these panels, we take p z = 0. 

In the present case, the p- wave component with S z = —1 is also induced (although we do not 
show this here), which is simply related to the S z = +1 component as §(p,S = 1 ,S Z = —1) = 
-$*(p,S = l,S* = l). 


n( /i'„) = Y 


i - /(w) - m P - q ,2) 


(39) 


p ^P+q /2 T ^p—q/2 W n 

Thus, the self-consistent solutions of the coupled equations (|26|) and (127)1 in Fig. |3](c) is valid 
for arbitrary values of the spin-orbit coupling strength A 2 . Since current experiments can 
reach the temperature region far below T c of a unitary Fermi gas in the absence of a spin-orbit 
interaction[l-4], the region in Eq. ()34l) is also accessible within the current experimental 
technique. In addition, since the superfluid order rapidly grows in the superfluid phase 
below T c , when the temperature can be lowered to some extent below T c , we would be 
able to obtain the p -wave condensate fraction, the value of which is close to that at T = 0 
obtained in this paper. 
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B. two-component spin-orbit interaction (Aj_ 7 ^ 0, A 2 = 0) 


When Aj_ 7 ^ 0 and A^ = 0, the spin-orbit interaction in Eq. (JTJ) consists of the cr x 
and a y components. In this two-component case, Eq. (TTB|) indicates that the p-wave pair 
amplitudes with S z = ±1 are induced. We show the detailed momentum dependence of the 
S z = 1 component in the upper and middle panels in Fig. [4j Although we do not show the 
S z = —1 component, it is simply related to the S z = 1 component as <h(p, 5 = 1,5', = —1) = 
—$*(p, 5 = 1, S z — 1). Thus, just after the s-wave interaction is suddenly replaced by the 
p-wave one in Eq. (I20jh one has the p-wave superfluid order parameters A p (p,S z = ±1). 
Because |A P (; p, S z — 1)| = | A p (p, S z = — 1)| , this p-wave superfluid phase is just the planar 


state 


41 


43]. Under the matrix representation in Eq. one has 


A(p) ~ 


-px + iPy 0 

0 Px + ip y 


(40) 


We briefly note that the planar state has not been realized in liquid 3 He. 

As in the single-component case, the pair amplitude dominantly appears around the 
Fermi surface in the weak-coupling BCS regime. When U s = 0, the present two-component 
spin-orbit interaction gives 




(p_i_ ± m\±)‘ 
2 m 


p 2 


(41) 


where /ij_ = p + mA^/2. At p, = 0, one obtains two Fermi surfaces, both of that are centered 
at p — 0, with the radii, 


k 


± 

F 


i/2m/ii ± mAi, p± > m\ 2 ± /2, 
+ mAi, p± < m\\/2. 


(42) 


Figures 0](al)-(cl) show that both the s-wave and p-wave pair amplitudes are large around 
these Fermi surfaces, as expected. 

In the two-component case, one cannot eliminate the A_L-dependence from the BCS- 
Leggett coupled equations (TTBT) and As a result, the s-wave superfluid order parameter 
A s , as well as the effective Fermi chemical potential jl ±, depend on the spin-orbit coupling 
strength A_l, as shown in Figs. Eta) and (c). Because of this, we see in Fig. Etb) and 
(d) that, not only the p-wave condensate fraction N C (S = 1,5~ = ±1), but also the total 
condensate fraction Ad is affected by the spin-orbit interaction, which is different from the 
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FIG. 5: (Color online) (a) Calculated s-wave superfluid order parameter A s and (c) effective 
chemical potential = n + mA|/2 at T = 0 in the case of two-component spin-orbit interaction 
(Ax 7 ^ 0, X z = 0). (b) p -wave condensate fraction N C (S = 1,S Z = ±1). (d) Total condensate 
fraction Ad. 


single-component case discussed in the previous subsection. The Ax dependence of Ad can 
also be confirmed by analytically carrying out the momentum summation in Eq. (1231) as 


N[ 

N 


37tA, 


16\/2i 


3/2 


, / 9 , A9 , VmX± vV + - mX 2 ± 2 

! l + V d + A 2 + —-— arccos - - 

2 Vdx + 


(43) 


Figure [5](d) shows that the total condensate fraction Ad in the BCS side ((&Fa s ) -1 ^ 0) 
is remarkably enhanced by the spin-orbit interaction, to be comparable to the value in the 
BCS regime with Ax = 0. Reflecting this enhancement, we find from Fig. ED(b) that the 
p-wave component N C (S = 1 ,S Z = ±1) is also enhanced in the BCS regime, when Ax is 
large. Thus, as shown in Fig. E)(a), the peak position of N C (S = 1 ,S Z = ±1) shifts to the 
BCS side, compared to the single component case shown in Fig. [31(b). 

Although the superfluid phase transition temperature T c is usually low in the weak- 
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(a) 



FIG. 6: (Color online) (a) Peak position of p -wave condensate fraction N C (S = 1 ,S Z = ±1) at 
T = 0 in the case of two-component spin-orbit interaction. The background intensity represents 
the magnitude of T c . For convenience, we also show the three-dimensional plot of T c in panel (b). 


coupling BCS regime, Fig. El)b) shows that T c around the “optimal region” (peak line in 
Fig. [6](a) ) in the BCS regime is also enhanced by the spin-orbit interaction. Since we need 
to reach the s-wave superfluid phase in order to prepare the p -wave pair amplitude, this 
enhancement of T c is favorable for our purpose. We briefly note that the enhancement of T c 
by a Rashba-type spin-orbit interaction has recently been pointed out 71]. 

The reason for the enhancement of the condensate fraction, as well as T c , in the weak- 


coupling BCS regime is the formation of two-body bound molecules[65|, IZlt 17314761]. which 


are also referred to as rashbons in the literature. In the present case, the degeneracy of the 
lowest energy level in the lower band in Eq. (PXTT) leads to a two-dimensional-like density 
of states around the bottom of this band. This “low-dimensional effect” stabilizes a two- 
body bound state (rashbon), even in the weak-coupling BCS regime 6hl. IfiI 1z3-0], where 
such a two-body bound state does not appear in the ordinary three-dimensional system. 
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FIG. 7: (Color online) Calculated effective chemical potential jl± = /r + mA^/2 at T c as a function 
of the spin-orbit coupling strength A_l. E^/2 = E^/2 + mX\/2 1 where E b is the binding energy of 
a two-body bound state, determined from Eq. (1441) . 


Then, the superfluid phase transition in the strong spin-orbit coupling regime is dominated 
by the BEC of rashbons, giving a high T c even in the BCS regime. This is similar to the 
case of the ordinary strong-coupling BEC regime of an ultracold Fermi gas 


29 


39], where 


the superfluid phase transition is dominated by tightly bound molecules that are formed 
by a strong pairing interaction. Indeed, as shown in Fig. [71 when the spin-orbit coupling 
strength A_l increases, the effective Fermi chemical potential fl± = /i + mX\/2 at T c becomes 
negative, to approach E h /2 = E h /2 + mX 2 L /2, where Ep is the binding energy of a two-body 


bound state, determined from the equation[65), [71, 


73 


76], 


1 = 


47ra, 


m 


E 


^E 


1 

2 £r 


(E h < -m\ 2 j2). 


(44) 


2 ^ | 2(^Ep T ciA_i_p_l) E\) i.cp 
Since the chemical potential physically describes energy to add a particle to the system, Fig. 
[7] indicates that most Fermi atoms form two-body bound molecules in the strong spin-orbit 
coupling regime. 
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FIG. 8: (Color online) Calculated p -wave condensate fraction at T = 0 in the case of three compo¬ 
nent spin-orbit interaction (Aj_, \ z ) = A(sin6h, cos 6\). We set (fcFa s ) _1 = 0. (a) N C (S = 1, S z = 0). 
(b) N C (S = 1 ,S Z = ±1). In this figure, 9\ = 0 and 9\ = 7t/2, respectively, correspond to the single¬ 
component and two-component spin-orbit interaction, discussed in the previous subsections. 


C. Three component spin-orbit interaction (Aj_ ^ 0, \ z ^ 0) 

When the spin-orbit interaction has all the a x , a y , and a z components [(Aj_, \ z ) = 
A(cos Ox, sin0 A ) (0 < 9\ < 7 t/2 )], all the p- wave condensate fractions, N C (S = 1, S z = ±1, 0), 
are induced, as shown in Fig. |8l 

When the spin-orbit interaction is isotropic (A_l = A 2 or 9\ = 7t/ 4), Figs. [91(a) and (b) 
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(a) 



FIG. 9: (Color online) (a) p- wave condensate fraction N C (S = 1 ,S Z = 0, ±1) at T = 0 in the case of 
isotropic spin-orbit interaction A_l = A^ = \/y/2 (0\ = tt/4). (b) Peak position of N C (S = 1 ,S Z = 
0, ±1) evaluated from panel (a). The background intensity shows T c , which is also shown in panel 
(c) for clarity. 


indicate that the BCS side is suitable for our purpose. Although we do not explicitly show 
the result here, the total condensate fraction, which is analytically given by 


C 

N 


3ttAo 


16^2 '< 


3/2 




i + 


mA 2 /4 

yP+Af 


(45) 
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(where fl = p + mA 2 /4) is also enhanced in the BCS side by the rashbon formation discussed 
in the previous subsection. In addition, T c is also enhanced by the same mechanism (See 
Fig. [21(c)), so that this regime is still considered to experimentally be accessible. In this case, 
when the s-wave pairing interaction is suddenly replaced by the p-wave one in Eq. (1201) . 
as the initial p -wave state, we can prepare the BW (Balian-Werthamer) superfluid order 


parameter 


41 


43], having the form 

A(p) ~ 


-Px + iPy Pz 
Pz Px + iPy 


(46) 


The BW state has been realized in superfluid 3 He 41]. 


IV. SUMMARY 

To summarize, we have discussed a possible idea to realize a p-wave superfluid Fermi 
gas. In contrast to the ordinary approach where one tries to cool down a Fermi gas with 
a strong p-wave pairing interaction, our idea consists of two stages. In the Erst stage, we 
only prepare a p-wave pair amplitude by, not using a p-wave interaction, but using the 
phenomenon that a p-wave pair amplitude is induced in a s-wave Fermi superfluid with 
broken inversion symmetry by an antisymmetric spin-orbit interaction. Then, in the second 
stage, we suddenly replace the s-wave interaction by a p-wave one, to produce the p-wave 
superfluid order parameter that is essentially given by the product of the p-wave interaction 
and the p-wave pair amplitude which has been prepared in the first stage. In this paper, we 
have assessed the first stage of our idea, by evaluating how large the p-wave pair amplitude 
can be induced in a spin-orbit coupled s-wave superfluid state. We clariEed the region where 
a large p-wave pair amplitude is obtained, in the phase diagram of a Fermi gas with respect 
to the strengths of the s-wave pairing interaction and the spin-orbit coupling. Within the 
framework of the NSR theory, we also confirmed that the s-wave superfluid phase in this 
optimal region is accessible within the current experimental technique. 

The key of our idea is that the pairing symmetry of a Fermi superfluid is just the symmetry 
of the superfluid order parameter, which is essentially given by the product of a pairing 
interaction and a pair amplitude. Thus, even when a p-wave pair amplitude is induced in 
a s-wave superfluid Fermi gas, the system is still in the s-wave superfluid state because of 
the vanishing p-wave interaction. This also implies the possibility that one may artificially 
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produce a p-wave Fermi superfluid by separately preparing a p-wave pair amplitude and a p- 
wave interaction. Our idea just uses this possibility, that is, we separately prepare these two 
quantity by using two sophisticated techniques developed in cold atom physics: the artificial 
gauge field and the tunable pairing interaction associated with a Feshbach resonance. 

At present, all the experiments aiming a p-wave superfluid Fermi gas uses a p-wave 
pairing interaction from the beginning. While this approach is straightforward, it seems 
suffering from the short lifetime of p-wave molecules, as well as the particle loss by the 
p-wave interaction. In contrast, since our idea does not rely on the p-wave interaction in the 
first stage, this difficulty can be avoided to some extent. In addition, since we can start from 
the situation with a finite value of the p-wave superfluid order parameter and a finite value 
of the p-wave condensate fraction, even when the p-wave interaction eventually destroys this 
superfluid state, we expect to be able to realize a p-wave superfluid Fermi gas for a while, 
after the p-wave interaction is introduced. 

In this paper, we have only discussed the first stage of our idea. Since the p-wave super- 
fluid state which is artificially produced would be in the non-equilibrium state, as the next 
step, we need to examine how this p-wave state relaxes into the thermodynamically stable 
state, after the p-wave interaction is turned on. In addition, since there are various p-wave 
superfluid phases, such as the BW phase and polar phase, the initial p-wave superfluid state 
may not be the most stable state for a given strength of a p-wave interaction and a given 
value of temperature. Thus, in order to soon reach the thermodynamically stable state, 
one need to clarify the optimal condition for the initial p-wave pairing state. In our future 
papers, we will discuss these interesting problems existing in the second stage of our idea. 
Since both the spin and orbital degrees of freedom are active in a p-wave Fermi superfluid, 
the realization of a p-wave superfluid Fermi gas would enable us to discuss much richer 
physics than the case of the s-wave superfluid Fermi gas. Since current experiments toward 
the realization of a unconventional superfluid Fermi gas faces various difficulties, our results 
would be useful for the exploration of a route to accomplish this exciting challenge in cold 
Fermi gas physics. 
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